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ABSTRACT 

The specific globular cluster frequencies (SV) for 24 dwarf elliptical (dE) 
galaxies in the Virgo and Fornax Clusters and the Leo Group imaged with the 
Hubble Space Telescope are presented. Combining all available data, we find 
that for nucleated dEs — which are spatially distributed like giant ellipticals 
in galaxy clusters — >S7v(dE,N) = 6.5 ± 1.2 and 5V increases with M v , while 
for non-nucleated dEs — which are distributed like late-type galaxies - 
^(dE, noN) = 3.1 ± 0.5 and there is little or no trend with My. The Sn values 
for dE galaxies are thus on average significantly higher than those for late-type 
galaxies, which have Sn ^ 1- This suggests that dE galaxies are more akin to 
giant Es than to late-type galaxies. If there are dormant or stripped irregulars 
hiding among the dE population, they are likely to be among the non-nucleated 
dEs. Furthermore, the similarities in the properties of the globular clusters and 
in the spatial distributions of dE,Ns and giant Es suggest that neither galaxy 
mass or galaxy metallicity is responsible for high values of Sn- Instead, most 
metal-poor GCs may have formed in dwarf-sized fragments that merged into 
larger galaxies. 

Subject headings: galaxies: star clusters — galaxies: elliptical and lenticular, cD 
— galaxies: nuclei 
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1. Introduction 

There are several distinct classes of galaxies with smooth surface brightness profiles 
and absolute magnitudes fainter than Mb = —18. Brighter than Mb = —15, most dwarf 
elliptical (dE) galaxies have central nuclei which are unresolved from the ground at the 
distance of the Virgo Cluster. However, there are quite a few dEs in this luminosity range 
that have no sign of a nucleus. Within the Virgo and Fornax clusters, the non-nucleated 
dEs form an extended population, with a spatial distribution similar to the spirals and 
irregulars rather than the bright ellipticals and the rest of the dEs (Ferguson & Sandage 
1989). A possible explanation for the difference is that the non- nucleated dEs are stripped 
or harassed dwarf irregular (dl) galaxies (e.g. Lin & Faber 1983; Moore et al. 1998). The 
most popular alternative scenario for dE formation (e.g. Dekel & Silk 1986) is that they 
formed in one monolithic collapse, and subsequently ejected their ISM via supernova winds. 
To distinguish these scenarios we require an observational probe of the earliest and most 
important episodes of star formation in dwarf galaxies. Globular cluster (GC) systems may 
provide such a probe. 

The specific globular cluster frequency (Sn, the number of clusters normalized to a 
galaxy with My = —15) is correlated with Hubble type such that it increases from late 
type spirals to early type spirals and ellipticals, suggesting a correlation with the size of the 
bulge component (Harris 1991a). The empirical variation of Sn from early to late types 
suggests a simple test for whether dEs are more closely related to giant ellipticals (Es) or 
to irregulars (Harris 1991b). If part of the E family, they should have Sn ~ 2 — 5; if 
part of the late-type disk galaxy family, we expect Sn ^ 1- Finally, one might expect that 
nucleated dEs will have higher values of Sn than non-nucleated dEs. 

We have begun a Hubble Space Telescope snapshot survey in order to study the 
properties of the GCs and nuclei of dE galaxies. In this paper we present results on Sn for 
24 galaxies in the Virgo and Fornax Clusters and the Leo Group. Table 1 gives the galaxy 
names morphological types from photographic catalogs. This sample represents all the data 
suitable for determining Sn taken during HST Cycle 6. We assume distance moduli of 
31.2 for the Virgo Cluster, 31.4 for the Fornax Cluster, and 30.3 for the Leo Group. When 
presenting an uncertainty of an average quantity we quote the standard deviation of the 
mean. The observations and analysis are summarized in Section 2. We present the results 
on the specific frequencies in Section 3 and discuss some implications in Section 4. The 
data will be more fully presented in Miller et al. (1998, Paper II). 
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2. Observations and Data Reduction 

WFPC2 snapshot images are taken in the F555W (2 x 230 sec) and F814W (300 sec) 
bandpasses with the galaxies centered in chip WF3. The procedures for object detection 
and classification, aperture photometry, and photometric calibration are similar to those 
used by Miller et al. (1997). From these data we can determine the (V — I) colors of GC 
candidates down to a completeness limit of V ps 25. 

The specific globular cluster frequency, Sn, is defined as Sn = N C 10 0A( - Mv+1 ^ where 
N c is the total number of clusters (Harris & van den Bergh 1981). Cluster candidates are 
selected based on color (0.5 < V — I < 1.5), size (FWHM < 2.5 pixels), and accuracy of 
photometry (err(V — I) < 0.3 mag). The net number of cluster candidates for a galaxy, N c , 
is the total number of cluster candidates in WF3 minus the expected number of background 
and foreground objects. The number density of background objects are determined 
from the objects in WF2 and WF4. We must also correct for undetected faint clusters. 
The luminosity function of all the cluster candidates is consistent with a Gaussian with 
M v eak —7.3 and a ~ 1.0. Tests then show that we are detecting 85-90% of the clusters in 
Virgo and Fornax (down to My ps —6) and ~ 98% in Leo (down to My ps —5). Therefore, 
we multiply the net number of detected GCs for Virgo and Fornax galaxies by 1.15 to 
correct for this incompleteness. 

We also measure the corresponding galaxy luminosity in WF3. After sky subtraction, 
most point sources and resolved background galaxies are removed from the image. The 
only cluster candidates which are not excised are the nuclei with My < —10.5 (see below). 
Then, the total counts remaining in the image are summed and corrected for foreground 
extinction to determine My. Uncertainties range between 0.1 mag for the brightest systems 
to ~ 0.4 mag for the faintest. 

An important problem is the classification of galaxies into nucleated and non-nucleated 
samples. Our first criterion for classifying nuclei is the projected radial distance, i? pro j, 
from the galaxy center as defined by the centroid of the light distribution after removing 
any central point-like sources. The uncertainties in the central coordinates are determined 
by measuring centroids with a range of centering-box widths. These uncertainties yield the 
number of standard deviations, a pro j, that the central objects are offset from the centers. 
As a further test, we have used Monte-Carlo simulations to compute the probability, P, 
that an object as bright as the central cluster candidate would be found within i? pro j of 
the center by chance. In each trial for a given galaxy, iV c artificial clusters are drawn from 
a Gaussian luminosity function and are given the same exponential radial distribution as 
measured for that galaxy. Many trials are run to reduce uncertainties due to small numbers. 
We classify a galaxy as nucleated if the center-most cluster candidate is within 3fjp r0 j of the 
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galaxy center and P ^ 1%. 

In most cases the nucleus turns out to be the most luminous GC, but in a few galaxies 
it is significantly brighter than the typical GC. With My <^ —10.5 these "clusters" are 
brighter than any globular cluster in the Galaxy or M31, and comparable to the brightest 
in M87 (Djorgovski 1993; Battistini et al. 1987; Whitmore et al. 1995). These objects may 
be more like true nuclei. It is unclear whether nuclei of large galaxies are kinematically and 
structurally distinct from star clusters or simply giant versions of the latter (Kormendy & 
McClure 1993; Carollo et al. 1997). To be conservative we have chosen not to count nuclei 
with My < —10.5 as clusters for the purposes of calculating Sn- 

3. Results 

Table 1 gives our measured values of My, N c (corrected for the faint end of the LF), 
Sn, -Rproj, Oproj by which -Rproj differs from the galaxy center, and probability P of finding 
the center-most cluster candidate within i? pro j by chance for 24 dE galaxies. The galaxies 
have been split into nucleated and non-nucleated categories based on the criteria mentioned 
above. In two cases our classification differs from the ground-based classification. FCC 324 
was originally in our non-nucleated sample but we find that it is nucleated. VCC 9 
has a catalog type of dElN. Its most luminous cluster is close to the center, but with 
.Rproj = 155 pc it is too offset to meet our criteria for being nucleated. Two other galaxies 
may be intermediate types. VCC 1577 and FCC 48 have their most luminous GC within 
3cr of their centers, but these galaxies do not fulfill the probability test. The central regions 
of VCC 9 and VCC 1577 are shown in Figure ID. 

We have one galaxy, VCC 1254, in common with the ground-based study of Durrell et 
al. (1996). Adjusting their value to our distance modulus yields Sn = 6.4 ± 3.0. Our result 
from Table 1 is Sn = 6.2 ± 2.2 and agrees very well with their measurement. Therefore, our 
method appears reliable. 

Figure |2| shows log(5 , Ar) versus My for all dE galaxies with measured Sn- In the 
figure, galaxies with Sn = have been assigned log (Sat) = 0.0. There is a clear trend of 
increasing Sn with decreasing galaxy luminosity for the nucleated galaxies (filled symbols). 
A weighted least-squares fit to the nucleated galaxies with My < —13.5 gives a slope of 
0.24 ± 0.04 (solid line). The non-nucleated galaxies (open and half-filled symbols) do not 
show such a trend, a fit these points gives a slope of 0.10 ± 0.05 (dashed line). The trend for 
the nucleated galaxies was seen by Durrell et al. (1996) but our larger sample now indicates 
that there is a difference in Sn between nucleated and non- nucleated dEs. 
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At magnitudes fainter than My = — 14 the large fluctuations in Sn caused by small 
number statistics may blur the distinction between nucleated and non-nucleated galaxies. 
An extension of the fit to the nucleated galaxies in Figure |2| passes through the points for 
both the Fornax and Sagittarius dwarfs, the dEs with the highest measured Sn- Even 
though the total extent of Sagittarius is uncertain because of tidal disruption by the Galaxy, 
M 54 is generally thought to be its nucleus (Sarajedini & Layden 1995). Fornax is usually 
considered non-nucleated, but it has a high Sn like the other dE,Ns and it has a metal-rich 
GC only ~150 pc from its center. The triangle labeled "LG dSph" represents the ~ 12 
non-nucleated Local Group dwarf spheroidals with Sn = 0. For now, it may be best to 
think of the dE,Ns as defining the upper envelope to the trend of increasing Sn with My. 

4. Discussion 

High resolution HST images of a large sample of dE galaxies provide an improved 
estimate of their specific globular cluster frequency, and an indication that Sn depends 
on galaxy type and luminosity. We confirm earlier ground-based estimates of a high Sn 
for dE galaxies (Durrell et al. 1996). Adjusted to our distances their sample of 9 dE,N 
and 3 dE,noN galaxies gives a unweighted mean of Sn = 5.1 ± 1.0. For comparison, our 
sample gives 5.3 ± 1.1. These values are consistent, and are significantly higher than the 
mean for late type galaxies Sn = 0.5 ± 0.2 (Harris 1991a). In addition, we find for our 
sample that 5jv(dE,N) = 7.5 ± 1.8 and Sjv(dE,noN) = 2.8 ± 0.7. The difference in S N 
between dE,N and dE,noN in Durrell et a/.'s sample is much less apparent, probably due 
to smaller sample size. The combined sample with My < —13.5 shown in Figure |2| yields 
/Sjv(dE, N) = 6.5 ±1.2 and 5jv(dE, noN) = 3.1 ± 0.5; the Sn for nucleated dEs is a factor of 
two higher than for non-nucleated dEs. 

This result is a further argument against the formation of dE galaxies as a class via 
the quiescent evolution of dl galaxies. The statistics of Sn in dls are not very good, but 
typical values appear to around Sn = 0.5 (Harris 1991a). Simple removal of gas from a 
dl galaxy would result in modest fading. For example, consider a galaxy that forms stars 
at a constant rate for 5 Gyr, whereupon it is stripped of gas and ceases forming stars. In 
the 5 Gyr after star formation ends the galaxy fades by ~1.5 mag in V. So, if it started 
with Sn = 0.5, it would end up with Sn ~ 2.0, which is at the low end of the observed 
range for the dEs. A longer period of fading would result in higher Sn, but the resulting 
surface brightness would tend to be fainter than measured in bright dE,Ns. Similarly, a 
model in which dEs in general represent the quiescent periods between stochastic bursts 
of star formation in dls (Tyson & Scalo 1988) appears inconsistent with the high average 
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value of S N measured for dE galaxies. The conclusion is that most bright dE galaxies, and 
the dE,Ns in particular, are not simply faded or dormant dl galaxies. 

The low mean Sn among non-nucleated dE galaxies suggests that any faded dls that 
exist in the cluster populations are likely part of this class. A number of the dE,noN 
galaxies in our sample have < 1, consistent with the values measured in dls. Also, our 
calculation above shows that faded dls would have Sn close to the mean Sn that we measure 
for the dE,noN class. A further argument in favor of this possibility is the observation 
that the non- nucleated dEs brighter than Mb = — 14 have an extended spatial distribution 
in the Virgo and Fornax Clusters like the spatial distribution of spirals and irregulars. It 
would be interesting to search for evidence of gas either in or around some of the brighter 
non-nucleated dE galaxies in these clusters. Also, stripped dls may still exhibit rotation. 

The colors and luminosity functions of the GCs in dEs are also similar to the old, 
metal-poor GC populations in giant ellipticals (Paper II). Therefore, GCs in both giant 
Es and dEs seem to have formed at similar times, with the same metallicities, and with 
the same efficiencies. It is not obvious why this should be so. The velocity dispersions of 
typical giant ellipticals are more than 250 km sec -1 , while for dE galaxies typical velocity 
dispersions are likely to be less than 100 km sec -1 (Bender & Nieto 1990). Thus it appears 
that galaxy velocity dispersion (and thus potential well depth) is not a key parameter. In 
addition, the metallicities of our dE galaxies, if adequately represented by their colors, are 
probably less than half solar, while the metallicities for many giant ellipticals exceed solar. 
Thus, the current metallicity of the galaxy also does not seen to play an important role in 
driving Sn- How can such similar GCs form in such different parent galaxies? 

This can be explained if most GCs are formed in dwarf galaxy sized "fragments" that 
latter merge into giant ellipticals (Searle & Zinn 1978; Harris & Pudritz 1994). The fact 
that E galaxies and dE,N galaxies inhabit the densest regions of clusters suggests that 
the cluster environment is somehow important for GC formation, for example through 
stronger interactions, higher external pressures, or higher star formation rates during the 
epoch of cluster formation (McLaughlin & Pudritz 1996; Elmegreen & Efremov 1997; 
Harris, Harris, & McLaughlin 1998). These environments would also be prone to the most 
merging. Additional evidence for this type of scenario is the accretion of the GCs from the 
Sagittarius dwarf into the halo of the Galaxy. Thus, accretion of GC-rich dwarfs may be a 
natural explanation of the metal-poor GCs in giant ellipticals and spirals. 

Turning to the trend with luminosity, our data suggests that, at least among nucleated 
dEs, less luminous galaxies have higher Sn- We have investigated several possible 
explanations for this phenomenon. First, winds may be more efficient at quenching star 
formation in the lower luminosity dEs. Second, many of the original clusters may have 
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merged into the bright nucleus. A 10 6 M Q cluster moving at 50 km sec -1 through an 
isothermal halo will decay from a radius of 3 kpc in less than a Hubble time (Binney & 
Tremaine 1987). We have calculated that Sn could decrease by a factor of about 1.5 if the 
nuclei brighter than My = —10.5 are merged GCs. However, these nuclei do not appear to 
be simple conglomerates of GCs. They are on average 0.06 ± 0.02 mag redder in (V — I) 
than the mean color of the GCs in those galaxies. If this is due to metallicity (Couture et 
al. 1990), then the nuclei are about 0.3 dex more metal rich than the GCs. Thus, additional 
nuclear star formation, perhaps as a result of interactions, may have occurred (also see 
Durrell et al. 1996). 

In summary, we have confirmed earlier estimates of a high globular cluster specific 
frequency in dE galaxies in nearby clusters of galaxies. Further, we have shown evidence 
for a difference in Sn between nucleated and non-nucleated dEs, and for a trend in Sn 
with luminosity (at least for the dE,N galaxies). The results reinforce the status of dE,N 
galaxies as the likely low-luminosity tail of the giant-E galaxy sequence, and suggest that 
the processes that formed both types of galaxies were related and probably influenced 
by the cluster environment. The present trends with luminosity and type are tentative 
and must be confirmed with larger samples. Fortunately, these are relatively inexpensive 
observations to make with HST, and we are engaged in a extension of the survey that will 
observe up to 30 more faint dEs and up to 30 dls to improve the statistical comparisons. 

We would like to thank Bill Harris, Frangois Schweizer, and Tim de Zeeuw for their 
comments and Jen Mack for help with the surface photometry. This research has made use 
of NASA's Astrophysics Data System Abstract Service. Support for this work was provided 
by NASA through grant number GO-6352 from the Space Telescope Science Institute, 
which is operated by AURA, Inc., under NASA contract NAS5-26555. 
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Table I. S N for Dwarf Elliptical Galaxies from HST 
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0.4 
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9.8 ±3.8 
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267.9 
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3.4 ±3.7 
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6.0 
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5.1 ±4.4 
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1.9 
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(^Corrected for undetected portion of the LF. 
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Fig. 1. — The central 35" of two dE galaxies with possible off-center nuclei. Each panel 
shows a physical region about 3 kpc across. 
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Fig. 2. — Trends of log^^r) with My for all dE galaxies with measured <SV- Circles are from 
the current work, triangles are from Durrell et al. (1996) after adjusting to our assumed 
distances, and squares are from the Local Group. The solid line gives the weighted fit to the 
nucleated galaxies (filled symbols) and has a slope 0.24 ±0.04. The dashed line gives the fit 
to the galaxies with offset nuclei (FCC 48 and VCC 1577) and no nuclei (open symbols) and 



